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ABSTRACT 

We have obtained single-epoch optical photometry for 201 quasars, taken from the FIRST Bright 
Quasar Survey, which span a wide range in radio loudness. Comparison with the magnitudes of these 
objects on the POSS-I plates provides by far the largest sample of long-term variability amplitudes for 
radio-selected quasars yet produced. We find the quasars to be more variable in the blue than in the 
red band, consistent with work on optically selected samples. The previously noted trend of decreasing 
variability with increasing optical luminosity applies only to radio-quiet objects. Furthermore, we do 
not confirm a rise in variability amplitude with redshift, nor do we see any dependence on radio flux 
or luminosity. The variability over a radio-optical flux ratio range spanning a factor of 60,000 from 
radio-quiet to extreme radio-loud objects is largely constant, although there is a suggestion of greater 
variability in the extreme radio-loud objects. We demonstrate the importance of Malmquist bias in 
variability studies, and develop a procedure to correct for the bias in order to reveal the underlying 
variability properties of the sample. 

Subject headings: quasars:general — methods:statistical 



surveys 

ing from long-term variability were ~ 0.25 magnitudes. A 
conference devoted exclusively to AGN variability - on all 
timescales and throughout the electromagnetic spectrum 
- summarized the status of the subject at the beginning of 
the last decade (Miller and Wiita 1991). Since that time, 
over a dozen major studies of optical variability have ap- 
peared (Table 1), covering timespans up to thirty years 
and including over 1000 objects. Most, however, include 
only optically selected quasars; Smith et al. (1993), Netzer 
et al. (1996), Sirola et al. (1998), and Garcia et al. (1999) 
are the exceptions with a combined total of 160 radio-loud 
quasars in their samples. Considerable disagreement re- 
mains regarding the correlation of variability amplitude 
with quasar power, redshift, and other quasar properties 
(see the references in Table 1). 

We are in the process of compiling a very large sample of 
bright quasars using the Very Large Array (VLA) 1 FIRST 
survey to select candidates from the POSS-I plates for op- 
tical spectroscopic followup. As part of this program, we 
have obtained current photometric B and R magnitudes 
for over 200 confirmed quasars, providing variability es- 

1 The VLA is part of the National Radio Astronomy Observatory which is operated by Associated Universities, Inc., under cooperative 
agreement with the National Science Foundation. 
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1. INTRODUCTION 

Variability on timescales ranging from hours to decades 
has been one of the defining characteristics of active galac- 
tic nuclei (AGN) since their discovery nearly forty years 
ago. These changes in the source flux received at Earth can 
arise from many causes: variability in the source's central 
engine, changes in a relativistic beam's orientation or ve- 
locity, changes in absorption along the line of sight to the 
active nucleus, gravitational microlensing, and interplane- 
tary or interstellar scintillation. As such, AGN variability 
provides an important diagnostic in studies of the physics 
of the central engine, the nuclear environment, the proper- 
ties of the material along the line of sight, and, ultimately, 
of the evolution of the AGN population. 

Twenty-five years ago, Grandi and Tifft (1974) pub- 
lished the first compilation of quasar variability studies. 
Citing 92 references to previous work, they listed photo- 
graphic and photoelectric photometry for a heterogeneous 
sample of 86 bright (B < 18) quasars, concluding that 
typical uncertainties in quasar optical magnitudes result- 
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Table 1 

Summary of Long-term Quasar Optical Variability Studies 
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timates on proper timescales of 10-45 years for a large 
quasar sample spanning a broad range in radio loudness. 
In section 2, we briefly describe the FIRST quasar survey 
and the selection from it of our variability sample. We 
go on to describe the photometric observations and their 
reduction (§3) as well as the procedures used to calibrate 
the photometry of the archival POSS-I data (§4). Section 
5 presents our results including a comparison with previ- 
ous work and an analysis of the dependence of variability 
amplitudes on radio loudness, redshift and luminosity, as 
well as the important effect Malmquist bias has on the dis- 
tribution of the magnitude variations. We conclude with 
a summary of our results. 



Observations were conducted through standard B and R 
filters; typical integration times were 100-600 s divided 
into two frames to facilitate cosmic ray rejection. Stan- 
dard stars from the list of Landolt (1992) were observed 
at least twice per night to solve for the extinction cor- 
rection. Conditions during the observations ranged from 
photometric to light cirrus. 

The data were reduced using standard IRAF routines 
for flat-fielding, cosmic ray removal, calibration and pho- 
tometry. In particular, the IRAF task PHOT was used to 
perform aperture photometry on both the quasar target 
and all other stellar objects of similar magnitude in the 
field. 



2. THE FIRST BRIGHT QUASAR SURVEY 

The FIRST Bright Quasar Survey (FBQS) has recently 
been described in detail by White et al. (2000). Briefly, 
the survey is based on a comparison of the radio catalog 
from the VLA FIRST survey (Becker et al. 1995; White 
et al. 1997) and the Cambridge Automated Plate Machine 
scans of the Palomar Observatory Sky Survey-I (POSS-I) 
plates (APM - McMahon and Irwin 1992; McMahon et al. 
2000). All optical objects brighter than E (red) magnitude 
17.8 (after corrections for extinction) 2 classified as stellar 
on at least one of the two POSS-I plates, which are bluer 
than O — E = 2.0 and lie within 1.2" of a radio source, are 
selected as quasar candidates. Spectroscopic followup has 
demonstrated that over 55% of these objects are quasars, 
with another ~10% classified as BL Lac objects; the re- 
maining sources are a mix of narrow-line AGN, HII galax- 
ies, and radio galaxies lacking emission lines. Over 700 
quasars and BL Lacs have been found in the first ^2700 
deg 2 of the survey; ~ 75% are newly discovered objects. 

For our followup program of optical photometry, 201 
quasars were selected at random from both the portion 
of the north Galactic cap survey reported in White et al. 
(2000) and its extension to the FIRST equatorial strip in 
the south Galactic cap (Becker et al. 2001). Suitability 
for the Right Ascension range of each scheduled observ- 
ing run was the only selection criterion; the magnitude, 
radio flux density, and redshift distributions of the chosen 
quasars are statistically identical to these distributions for 
the FBQS sample as a whole. Flux densities for the sample 
cover a range 1.0 - 15,000 mJy, and objects with redshifts 
from 0.1 to 3.4 are included. The ratio of 5 GHz to 2500 A 
flux densities R* used to characterize the radio loudness 
of an object (Stocke et al. 1991) spans the range 0.3 to 
22,000; objects with R* < 10, generally characterized as 
radio-quiet quasars, constitute 46% of the sample. 

3. OPTICAL PHOTOMETRY 

Observations of the selected objects were carried out by 
one of us (RPSS) in a number of observing runs at the Lick 
1-m telescope over the period 20 December 1995 through 
6 June 1997. The 2048 x 2048 pixel Orbit CCD provides 
a 6' field of view, sufficient to include a number of stars 
in each image for use in calibrating the APM magnitudes. 

2 The initial selection for the FBQS came from the APM magnitudes. Subsequent calibration of these magnitudes using the APS scans (see 
McMahon et al. 2000 for details) led to a small number of objects already spectroscopically confirmed as quasars being dropped from the 
sample presented in White et al. (2000) because their revised POSS-I magnitudes fell slightly below the R = 17.8 threshold. This program, 
however, was begun before the recalibration. Thus twenty-one of the sources observed have POSS-I magnitudes slightly fainter than R = 17.8 
as can be seen in column 9 in Table 2. Eight of these have been published in Gregg et al. (1996), while the remainder appear here for the first 
time. 



4. CALIBRATION OF THE APM MAGNITUDES 

The preliminary photometric calibration of the APM 
POSS-I scans currently in use is described in some de- 
tail in McMahon et al. (2000). For stellar objects, the 
global rms photometric uncertainty is 0.5 magnitudes in 
the range of interest here (15 < O, E < 19), a precision in- 
sufficient for detecting the expected variability amplitudes 
for most quasars. Thus, we began by using our CCD pho- 
tometry in the 201 quasar fields to calibrate each of the 53 
POSS-I plates on which our objects were found. 

There are from 1 to 15 quasar fields per plate. In each 
field, we selected isolated stellar objects with 13 < B, R < 
19 for use in forming a sample for photometric calibra- 
tion. The magnitudes in the two CCD exposures for each 
star were compared; if they differed by more than 0.2 
magnitudes, the star was discarded. The stars were then 
matched to the APM catalog with a matching criterion 
of 10" (in order to account for proper motion over the 
~ 45-year interval between the POSS-I and our observa- 
tions); any star with no match was discarded. If two stars 
matched within this radius, the nearer one was chosen un- 
less it differed in brightness by more than 2 magnitudes. 
These procedures yielded a database of 1695 stars. 

We then performed an interative, linear least-squares fit 
to the magnitude differences (B — O, R—E) as a function of 
O, E magnitude in order to derive a magnitude-dependent 
calibration for each pair of plates; stars lying more than 
3(7 from the best fit on either plate were deleted from the 
database in order to eliminate variable stars, mismatched 
stars, images affected by cosmic ray hits, etc. We required 
a minimum of five stars per plate for an acceptable cali- 
bration; the median was 16 stars per plate, with the best 
covered regions having more than 80 stars per plate. The 
mean (and median) rms values from the fits to both sets of 
plates are 0.15 magnitudes; the fit for the POSS-I O plate 
1146 is displayed in Figure 1. The corrections were then 
applied to the O and E magnitudes for each quasar. We 
estimate that the combined uncertainty in the magnitude 
differences due to measurement and calibration errors is 
0.17 magnitudes for both the red and blue. 

As an additional check on the calibration procedure, we 
compared our CCD-corrected magnitudes with the APS- 
corrected magnitudes which we have derived for the entire 
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Fig. 1. — The calibration results for the POSST O plate 1146. Stars selected from our CCD images for use in the calibration arc represented 
by dots, while quasars are shown as asterisks. The rms for the 78 stars is 0.14 mag; note the substantially higher scatter for the quasars, as 
well as the marked tendency for the quasars to be fainter (positive B — O) at the current epoch. 



APM/FffiSTidentification catalog described in McMahon 
et al. (2000). The distribution of corrections was very sim- 
ilar; since the scatter in the APS magnitudes was slightly 
larger, we used our CCD plate-by-plate calibrations when 
they were available. For eight of our quasars, there were 
too few calibration stars per plate, and the APS magni- 
tudes were used; these sources are flagged in Table 2. 

5. LONG-TERM QUASAR VARIABILITY 

5.1. Photometric results 

In Table 2, we present the calibrated O, B, R, and E 
magnitudes for our 201 quasars. Each source's FIRST co- 
ordinates are given in the first two columns. The next 
two columns list the dates of observation for the POSS I 
plates (both colors were taken on the same night) and for 
our CCD observations, followed by the calibrated magni- 
tudes and B — O, R — E variations (positive B — O, R — E 
values mean the source is fainter now than when it was 
first detected on the POSS I plates). We also include the 
FIRST peak and integrated 20cm flux densities, the mea- 
sured redshift (from White et al. (2000), Becker et al. 
(2001), Gregg et al. (1996) and our own unpublished val- 
ues for the remaining objects with E > 17.8), and val- 
ues for the absolute blue magnitude, radio luminosity, and 
logi?*, the radio loudness parameter; an H = 50 km s _1 
Mpc -1 , O = 1, A = cosmology is assumed throughout. 
Note that the radio luminosities are for the radio core only 
and do not include flux from any extended lobes. The fi- 
nal column contains a radio morphology flag; about 10% of 



the sources are classical double-lobed sources, while most 
of the remainder are point-like at the ~ 5" resolution of 
the FIRST survey. 

By far the most variable object in the sample is the 
EGRET source 4C38.41 (FIRST J131059.4+323334), a 
well-known Optically Violent Variable which was 4.45 and 
4.17 magnitudes brighter in the E and O bands, respec- 
tively, in 1950 than when we observed it in 1996. Since 
Barbieri et al. (1977) saw this object change from B=15.85 
to B=19 over the course of four years, the large varia- 
tion we report is not unexpected. The only other object 
which varied by more than 1.5 magnitudes in cither band is 
the low-redshift (z — 0.28), core-dominated, flat-spectrum 
quasar B1719+357 for which B-0= 1.98. 

In all, 30 objects (15%) varied by at least a factor of 
2 in at least one band. Of these, only four have been 
reported previously in the literature as optical variables. 
For the sample as a whole, there is a large bias toward 
objects having faded since their POSS-I detections; a de- 
tailed analysis of this effect, a manifestation of Malmquist 
bias, is given below, following which we present the re- 
sults of comparing long-term variability with other quasar 
properties. 

5.2. Skewness of the 8m distribution 

Figure 2 shows the distributions of the change in quasar 
magnitude (dm = mcco — rnposs) m the red and blue 
bands. The distributions are clearly skewed, with many 
more quasars found to be dimmer in the CCD data than 
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Fig. 2. — Observed distribution of magnitude change 8m = mcCD — mposs f°r quasars in our sample. Objects on the right were fainter 
in the current CCD observations. The histogram is skewed due to Malmquist bias (sec text for discussion.) The dots are the observed counts, 
and the solid line shows the predicted counts in bins (with la uncertainty error bars) from a symmetrical, Malmquist-corrected model. The 
dashed line is a Gaussian distribution having the same rms as the model distribution, (a) Red magnitude differences (riiR — ms). The 
Gaussian width is a = 0.33. Error bars account for the use of these same data to compute the model values and show the expected range for 
data points about that mean, (b) Blue magnitude differences (mg — mo). The model (both histogram and Gaussian) arc derived from the 
red delta-magnitude distribution using the relation B — O = 1.3(R — E). In this case the model is derived completely independently from the 
data, so the error bars show the expected range for these data about the mean due to Poisson statistics. 



6 



Long-Term Quasar Variability 



are found to be brighter. This asymmetry is the result of 
Malmquist bias (Malmquist 1924) in the FBQS sample. 
Because the quasar number counts increase steeply with 
magnitude, a disproportionate fraction of quasars near the 
FBQS magnitude limit are included in the sample because 
small brightness fluctuations at the time of the POSS-I 
observation moved them above the optical threshold for 
inclusion in the sample. 

It is relatively straightforward to correct the observed 
dm distribution, P i ]s {8m) 1 for Malmquist bias and thus 
determine the true distribution, P true (5m). Let m be the 
true (current-epoch CCD) magnitude of the quasar, m 
be the magnitude brightness limit of the FBQS sample, 
and N(m) be the cumulative number of quasars brighter 
than m. For 5m > 0, the quasars were brighter on the 
POSS-I plates than they are currently, so the FBQS sam- 
ple was effectively looking deeper in the N(m) distribution 
and thus included more objects. The observed magnitude 
distribution is 



P bs(Sm) = 



Ptrue(5m)N(m + 5m) 
I-oo Ptrue(5m)N(m + 5m)d5m 



(1) 



The factor in the denominator simply normalizes P b s so 
that its integral is unity. 

For quasars at V ~ 17.5, the cumulative number distri- 
bution is approximately exponential: 



N(m) = jVoe( m - m °) lnF 



(2) 



where F is the factor by which the counts increase for a 
1 magnitude change (F ~ 7.6, La Franca and Cristiani 
1997). Thus the observed distribution is 



P obs (8m) = 



Ptme{5m)e 



5m In F 



JZ Ptrue(5m)e Sm ^ F d5m 



(3) 



This integral can be evaluated analytically for a Gaussian 
distribution. If P trU e — exp(— 8m 2 /2a 2 ) then 



Pobs{5m) = exp(-[5m - a 2 lnF] 2 /2cr 2 ) 



(4) 



A Gaussian distribution simply has its mean shifted by 
cr 2 lnF = 2.03cr 2 for F = 7.6. 

More generally, it is safe to assume that the true distri- 
bution of 8m is symmetrical about zero. Using this (weak) 
assumption, we can symmetrize the observed histogram of 
5m. Consider two histogram bins symmetrically placed 
about zero: \5m — 5mh\ < b and \5m + 5mh\ < b, where 
the bins are centered at ±5mh and the bin width is 2b. Let 
the number of objects in each 5m bin be N+ and N-. By 
the symmetry assumption, the true distribution has the 
same expected number of objects in each bin, N. If the 
bin width is small, the expected number of objects E(N) 
in the observed distribution is modified by a Malmquist 
bias factor W: 

E(N+) = NW+ = Nexp(5m h InF) , and (5) 
E(N_) =NW- = Nexp(-8m h \nF) . (6) 

The most accurate estimator of the true number of objects 
in the bin is then 



N = 



N + + N_ 
W+ + W- 



Some straightforward algebra will convince the reader that 
this gives the correct answer in the case of a Gaussian 
distribution (except for a constant normalization factor, 
which can be computed _directly from N.) The uncer- 
tainty in the estimator N is easily determined, since it 
has a Poisson noise distribution scaled by the W+ + W- 
factor. 

A subtle point is that the observed blue magnitude dif- 
ferences are biased only through their coupling to the red 
magnitudes, because the red magnitude was used to de- 
fine the sample. There is a weak bias that directly af- 
fects the blue due to a limit O — E < 2 on the colors of 
FBQS quasars. Thus, quasars that are near the FBQS 
limit E = 17.8 must be brighter than O = 19.8 to be 
included in the sample. However, most quasars are much 
bluer than this limit and so the bias in the blue 8m b distri- 
bution (which is apparent in Fig. 2) appears only because 
the blue variation is highly correlated with the red varia- 
tion. 

To understand how the bias affects the blue 8m b distri- 
bution, consider a hypothetical population of objects that 
do not vary at all in the red but do vary in the blue. For 
simplicity we also assume measurement errors are zero. At 
the FBQS discovery epoch, we simply select all the objects 
brighter than the sample magnitude limit. Some of those 
objects are brighter than average in the blue and some are 
fainter, but all are included in the sample because there 
is no discrimination against faint blue magnitudes. At the 
current epoch we would find that the blue magnitudes dif- 
fer but the red magnitudes are unchanged. The observed 
blue 8mB distribution would consquently be symmetrical 
about zero, with no Malmquist bias at all. 

Now consider another extreme hypothetical population, 
where both the red and blue magnitudes of the objects 
vary but the blue magnitude is completely determined by 
the red magnitude, such that 5ms = a5mR where a is 
a constant. In this case, the selected sample is biased 
by brightness fluctuations that bring objects with brighter 
R magnitudes into the sample; those objects also have 
brighter B magnitudes. At the current epoch, we there- 
fore find that many objects have dimmed in both B and R 
according to the hypothesized linear relationship. In this 
case the observed 5mn distribution is related to the true 
distribution as derived in equations (l)-(7) above, and the 
5mB distribution is merely a rescaled version of the 8mR 
distribution. For example, if a = 2, then the 5ms distri- 
bution will be twice as wide, shifted twice as far from zero, 
and (to conserve counts) half as high as the 5mu distribu- 
tion. This result is not the same as if the blue magnitudes 
were themselves Malmquist-biased; in that case, when the 
distribution is twice as wide, it is shifted four times as far 
from zero (Eq. 4) rather than twice as far. 

The real case is, of course, between these two ex- 
tremes, with the blue variation having components that 
are both correlated with and independent of the red vari- 
ation. Fig. 3 shows that 5ms is in fact strongly corre- 
lated with 5mR. The best-fit line through the points, 
allowing for measurement errors in both coordinates, is 
5ms = 1.30 5mR. If we assume that the uncertainty in 
SmR = 0.17 (an estimate of the combined measurement 
errors in E and R), then the scatter about the best fit line 
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Fig. 3. — A comparison of the R — E and B — O variations for our quasar sample. The dashed line indicates identical variations in the two 
colors, while the solid line has the best fit slope of 1.30. The most highly variable source (Siur ~ 4.5, Stub ~ 4.2) falls outside the plot. The 
magnitudes arc strongly correlated, with higher variability amplitudes in the blue. 



indicates that the uncertainty in 5ms — 0.24 in order that 
the x -sc t uare °f the fit be approximately 1 per degree of 
freedom. Since the measurement errors on B and O are 
similar to those in the red, the implied "intrinsic" scatter 
in Stub — 0.17. In other words, the portion of the blue 
variability that is uncorrelated with the red variability is 
relatively modest and is of the same order as the measure- 
ment error. 

Figure 4 shows the 5m distributions for both the red and 
blue symmetrized and corrected for the Malmquist bias 
(see Eq. 7). The predictions from the symmetrical model 
are in good agreement with the data. Note in particular 
the reasonably good agreement with the blue distribution; 
since the model in this case was derived from the red Sttir 
data, there are no free parameters in the fit. The red dis- 
tribution was converted to a blue model simply by scaling 
the 8m values by a factor 1.30, as indicated by the fit in 
Fig. 3. A Gaussian with the same rms as the model dis- 
tribution is also shown and appears to be quite a good 
representation for the bulk of the distribution, although 
there is certainly a non-Gaussian tail of high-amplitude 
variables. 

The predictions of the symmetrical model (and the as- 
sociated Gaussian) are also shown in Fig. 2 where they are 
overlain on the observed (biased) counts. This predicted 
distribution is obtained by multiplying the true symmetri- 
cal distribution by the Malmquist factor W from Eq. (4). 
As was mentioned above, the biased Gaussian distribu- 
tion is just a shifted Gaussian with the same width. Note 



that the Gaussian in Fig. 2a has only one free parameter 
(the width) because the shift is specified by the width and 
the normalization is determined by the total number of 
quasars; thus the good agreement with the observational 
distribution supports both the general approach adopted 
here and the idea that there is a roughly Gaussian under- 
lying distribution of amplitude fluctuations. 

In the following section we examine the dependence of 
variability on various quasar properties (e.g., redshift and 
radio- loudness). Malmquist bias affects all of these mea- 
surements. The correction procedure derived above can 
also be applied to any subset of the data to determine a 
Malmquist-corrected rms value for the variability of that 
subset. The approach we have adopted is to compute the 
symmetrized distribution for the subset (typically objects 
having some property falling within a moderately narrow 
range) and then calculate the rms from that distribution. 
From the analysis of the whole sample above, this ap- 
proach is seen to give a good measurement of the width of 
the (approximately Gaussian) histogram. We exclude the 
single object that varied by 4 magnitudes (4C38.41) from 
these calculations since it produces a wild outlying value 
for the rms of any bin into which it falls. 

5.3. Correlation with other quasar properties 

Figure 3 compares the variability amplitudes for the red 
and blue bands. The correlation between the two bands is 
apparent, as is the fact that the average variability ampli- 
tude is higher in the blue. This trend has been seen in pre- 
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Fig. 4. — The distribution of magnitude differences corrected for Malmquist bias under the assumption that the true distribution is sym- 
metrical about zero. As in Fig. 2, the dots are the observed counts corrected for the Malmquist factor in each bin, the solid line is the true 
distribution of Sm, and the dashed line is a model Gaussian; in both panels, the models are derived from the red data, (a) Red magnitude 
differences, (b) Blue magnitude differences. 
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vious optically-selected samples; e.g., Figure 2 in Trevese 
et al. (1997) suggests that for a factor of two change in 
wavelength, Stub — Sinn ~ 0.11. For our sample, the mean 
difference between the variability amplitudes is 0.13 mag. 
When the Malmquist bias is removed, the rms variations 
in the blue and red are 0.45 mag and 0.33 mag, respec- 
tively (Fig. 2). If we subtract in quadrature the 0.17 mag 
rms contributed by calibration errors (§4), the quasar blue 
and red rms values are 0.42 and 0.28 magnitudes. Thus the 
bias-corrected rms difference is 0.14 mag, in good agree- 
ment with Trevese et al. 

In Figure 5, we show the distribution of 8m as a func- 
tion of rcdshift. We have made no correction for emis- 
sion lines in our analysis; in that B, O and R, E are not 
identical bands, strong lines entering the bands at slightly 
different redshifts could be a cause for concern. We have 
marked the locations at which Mgll, CIV, and Lya enter 
and leave the POSS-I bands; no significant anomalies are 
seen. While it is possible that a few objects with very 
strong lines could have their 5m values affected by line 
emission, it does not appear that the statistical conclu- 
sions derived from our sample will be biased by this effect. 
Note again the strong preference for positive 5m values 
induced by the Malmquist bias. However, the amplitude 
of the variability in both bands is independent of z. Bin- 
ning the objects in rcdshift intervals of 0.5, we find, for 
example, that < 5mn > differs by only 0.01 mag in the 
z < 0.5 and z > 3 bins, and in none of the seven bins is the 
mean value more than 1. 5a from the overall mean value 
of 0.33 when correction for the Malmquist bias is applied 
(see above). These results are inconsistent with the ma- 
jority of recent studies on optically selected samples (e.g., 
Fcrnandes, Aretxaga, and Terlevich 1996 and references 
therein). Given the complex interaction of k-corrections, 
wavelength-dependent variability, the redshift-luminosity 
correlation in a flux- limited sample, etc., it is difficult to 
assess the significance of this disagreement. 

The range of redshifts in the sample (and the range of 
time over which the data at the two epochs were collected) 
means the proper time intervals sampled range from 10 to 
45 years. The mean values of the variability amplitudes 
for proper times less than, and greater than, 20 years are 
identical to within 0.02 and 0.01 magnitudes, the respec- 
tive errors in the means for the blue and red bands. 

Figures 6a and b illustrate the size of the effect 
Malmquist bias can introduce when searching for correla- 
tions between variability amplitude and other quasar prop- 
erties. We show the median amplitudes over a span of 10 7 
in radio luminosity both before and after correcting for 
the Malmquist bias using the prescription outlined above. 
While before correction, the plots are noisy but suggest 
a positive correlation with luminosity in the R band, the 
corrected values are flat, with all points consistent with 
the median value to within less than 2a. We also find 
that variability amplitude is uncorrelated with optical ap- 
parent magnitude when Malmquist-corrected distributions 
are employed. 

The high resolution of the FIRST survey allows us to 
comment on the radio morphology dependence of the vari- 
ability. The FBQS selection criteria require a radio com- 
ponent coincident with the optical counterpart, so our sur- 

3 Comparing the radio quiet sources with Mb < —28 to those in 
at the 6.1cr level. 



vey may be partially incomplete for sources with very weak 
core components and dominant radio lobes (classical FR-II 
radio doubles), although our recent extension of the FBQS 
to include candidate doubles suggests the incompleteness 
is less than 5% (Becker et al. 2001). Furthermore, it is 
possible that the lobe components for some of our identi- 
fied quasars have been resolved out in our high-resolution 
radio observations, although, again, the fraction of such 
objects is small, since most have not been found to have 
extended components in the much lower-resolution NVSS 
survey (Condon et al. 1998). 

In total, roughly 10% of the sources in the current sam- 
ple show extended emission, and most of these are FR-II 
objects. We have computed the rms(i3 — O) and rms(i? — 
E) values for the single and complex sources separately 
and, surprisingly, find that those sources with extended 
compnents are slightly more variable; e.g., rms(i? — E) = 
0.320 ± 0.010 for the single sources, while rms(i? — E) = 
0.376 ± 0.028 for the complex objects. This result is not 
highly significant (~ 2a) but is in the opposite sense ex- 
pected from standard unified models of AGN in which the 
sources with extended lobes are oriented roughly normal 
to the line of sight, while sources with bright core compo- 
nents have their jets beamed toward us and should thus, on 
average, be more variable. Given possible incompleteness 
in the sample and the low significance of the distinction, 
it is clear that a larger sample of FR-II sources would be 
required to draw quantitative conclusions. 

Fig. 7 shows the dependence of variability amplitude 
on the radio-loudness parameter over a range in R* of 
10,000. The curves in the two bands are remarkably simi- 
lar, and suggest a nearly constant amplitude over most of 
the range, with a possible turn up for the extreme radio- 
loud objects. 

The final figure (Fig. 8) shows a three-dimensional dis- 
tribution of variability amplitude as a function of both 
R* and absolute magnitude. Fluctuations from the mean 
value of less than la have been suppressed (by shrinking 
all values toward the mean) to emphasize the more sig- 
nificant features of the distribution. The trend noted by 
Hook et al. (1994) in their optically selected sample for 
variability amplitude to decrease with increasing absolute 
magnitude is apparent here for the objects with R* < 10 
(the radio-quiet portion of our sample) 3 , but disappears 
for objects with higher radio-to-optical flux ratios. The 
dependence of this trend on radio loudness could arise if 
the overall variability is composed of two parts: one due 
to long-term changes in the structure of the accretion disk 
giving rise to changes in the thermal tail of the 'big blue 
bump', and the other arising from fluctuations in the non- 
thermal, beamed emission. The most luminous objects 
with large, slowly varying disks would have a smaller am- 
plitude of variability unless the nonthermal emission comes 
to dominate the optical flux; the higher amplitude variabil- 
ity for the radio-loudest objects, apparent in the figure, is 
consistent with assigning the greater variability amplitude 
to the beamed emission. 

The decrease in variability seen here for the lowest lu- 
minosity objects is significant at the 4.9<x level (comparing 
M B > -24 to -24 > M B > -28) and has not been re- 
range —28 < Mb < —24, the rms variability amplitudes are different 
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ported previously as a consequence of the exclusion from 
most quasar samples of objects with Mb > —23.5. Fur- 
thermore, the trend apparently does not continue toward 
lower luminosity Seyfert galaxies; in the compilation by 
Winkler et al. (1992), six of the seven highest amplitude 
variables over a four-year period are drawn from the up- 
per half of the luminosity distribution (median luminosity 
M B = -21.2). 

6. SUMMARY AND CONCLUSIONS 

We have presented a study characterizing the long-term 
(~ 45 yr) variability properties for a radio-selected sample 
of 201 quasars spanning large ranges in redshift, luminos- 
ity, and radio loudness. We find little or no dependence 
of the variability amplitude on radio or optical flux, radio 
luminosity, or redshift. We do see higher mean amplitudes 
in the blue band than the red band, consistent with the 
reported wavelength dependence in optical samples. The 
dependence on radio loudness is largely flat, with a slight 
upturn for the extremely radio-loud objects. The decrease 
in variability amplitude for optically luminous quasars is 
confined to radio-quiet objects. We also demonstrate the 
importance of considering Malmquist bias when charac- 
terizing the variability properties of a quasar sample, and 
describe a technique for correcting observed distributions 
of variability amplitude. 

The complex interdependence of the observed variability 
on frequency, redshift, bolomctric luminosity, radio loud- 



ness, and possibly other parameters makes it difficult to 
draw definitive conclusions regarding the physics of the 
central engine or quasar evolution from a study such as 
this. A larger study of a carefully selected sample of FBQS 
objects, supplemented by a sample of much fainter FIRST 
counterparts, could provide the breadth of coverage in lu- 
minosity over a broad range of redshifts necessary to begin 
to untangle these effects. 
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Fig. 8. — Red variability amplitude is displayed as a function of the radio loudness parameter R* and the absolute blue magnitude. The 
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Table 2 — Continued 
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17. 


.44 


0. 


02 


16. 


61 


16. 


81 


-0 


.17 


46 


.17 


48 


.36 


1 


.741 


-28. 


.1 


33 


.5 


1 


.79 


11 03 


13 


.31 


+30 


11 


42 


.9 


1955 


304 


1997. 


.189 




18 


18 


17. 


.65 





.53 


17. 


64 


17. 


90 


-0 


.26 


107. 


,68 


108 


.99 





.387 


-24. 


3 


32 


.6 


2 


.37 


11 06 


05 


.69 


+30 


51 


08 


.5 


1955 


.304 


1997. 


.189 




17 


.45 


17 


10 





35 


16. 


86 


16. 


80 


0. 


06 


1 


.69 


1 


.92 


1 


638 


-28. 


3 


32 


.0 





.26 


11 10 


40 


.21 


+30 


19 


09 


.9 


1955 


.304 


1997. 


.189 




18 


.97 


17. 


.98 





.99 


18. 


21 


17. 


86 





38 


20 


.08 


22 


.77 


1 


.522 


-27. 


.3 


33 


.1 


1 


.69 


11 22 


41 


.49 


+30 


35 


35 


.1 


1950. 


.425 


1997. 


.189 




17 


65 


17. 


.01 





64 


16. 


74 


16. 


74 


-0 


.00 


10 


.06 


10 


.02 


1 


.809 


-28. 


,7 


32 


.8 





.93 


11 31 


09 


.29 


+26 


32 


08 


.1 


1950. 


.425 


1997. 


595 




16. 


96 


16 


.41 


0. 


52 


16. 


40 


15. 


46 





.94 


2 


.00 


2 


.72 





.243 


-24. 


.5 


30 


.6 





31 


12 02 


43 


.54 


+37 


35 


52 


.0 


1956. 


.430 


1997. 


578 


* 


18 


48 


17 


.54 


0. 


94 


17. 


77 


17. 


36 





.11 


16 


.29 


17 


.24 


1 


.194 


-27. 


.1 


32 


.7 


1 


.43 


12 04 


37 


.56 


+28 


51 


25 


.1 


1955 


373 


1997. 


.578 




17 


36 


17 


.07 





.29 


16. 


64 


16. 


51 





13 


6 


.24 


5 


.57 


1 


.129 


-27. 


.5 


32 


.2 





.80 


13 10 


59 


.41 


+32 


33 


34 


.6 


1950. 


.449 


1996. 


.279 




20. 


96 


16 


.79 


1. 


17 


20. 


56 


16. 


11 


4. 


45 


243. 


.42 


245 


.97 


1 


.650 


-28. 


.7 


34 


.2 


2 


.24 


13 11 


31 


.98 


+37 


31 


08 


.1 


1950 


449 


1996. 


.304 




17 


.31 


17. 


.11 





.20 


17. 


11 


16. 


83 


0. 


28 


1 


.20 


1 


.53 


1 


.272 


-27. 


,7 


31 


.7 





.19 


13 12 


17 


.75 


+35 


15 


20 


.6 


1950 


449 


1996. 


.304 




15 


55 


15 


.66 


-0 


.11 


15 


05 


14. 


86 


0. 


19 


43 


.92 


44 


.72 





183 


-24. 


,6 


31 


.5 


1 


.22 


13 14 


23 


.14 


+29 


10 


01 


.7 


1955. 


367 


1996. 


304 




18. 


73 


18. 


.30 


0. 


43 


17. 


90 


17. 


.52 





.38 


7 


.53 


7 


.22 


1. 


.094 


-26. 


.2 


32 


.3 


1 


.38 


13 26 


29 


.06 


+36 


23 


31 


.3 


1950 


452 


1996. 


696 




18 


10 


17 


.60 





.50 


17. 


63 


17. 


23 


0. 


40 


1 


.75 


1 


.07 





389 


-24. 


,4 


30 


.8 





.56 


13 31 


08 


.31 


+30 


30 


32 


.9 


1950. 


.507 


1996. 


304 




17. 


49 


17. 


38 





11 


17 


11 


17. 


11 





03 


14777. 


.96 


15023. 


.95 





.848 


-26 


.5 


35 


.4 


4. 


.34 


13 33 


09 


.71 


+36 


20 


30 


.4 


1950. 


452 


1996. 


701 




17. 


95 


17 


.46 


0. 


49 


18. 


53 


17. 


38 


1 


.15 


16 


.63 


17 


.50 





.679 


-25. 


.9 


32 


3 


1 


.16 


13 34 


03 


.81 


+37 


01 


01 


.0 


1950. 


.452 


1996. 


.696 




18 


.38 


17. 


.65 





.73 


17. 


45 


17. 


05 





.40 


13 


.17 


13 


.31 


1 


.765 


-28. 





32 


.9 


1 


.31 


13 41 


59 


.89 


+37 


07 


10 


.5 


1950 


.452 


1996. 


701 




17 


16 


16 


.67 





.49 


17. 


84 


16. 


61 


1 


23 


24 


.50 


25 


.42 


1 


.106 


-27. 


,8 


32 


.8 


1 


.26 


13 42 


51 


.40 


+28 


28 


06 


.1 


1950. 


507 


1996. 


515 




18. 


50 


18. 


.20 


0. 


30 


17. 


83 


17. 


77 





06 


65 


.57 


68 


.41 


1 


.040 


-26. 


1 


33 


.2 


2 


.30 


13 43 


00 


.18 


+28 


14 


07 


.3 


1950. 


.507 


1996. 


515 




17. 


23 


17. 


.00 


0. 


23 


16 


77 


16 


43 





.34 


217. 


.37 


246 


.65 





.905 


-27. 


.0 


33 


.7 


2. 


.40 


13 48 


01 


.37 


+28 


40 


25 


.6 


1958. 


.373 


1996. 


515 




18 


00 


17. 


.89 


0. 


11 


17. 


77 


17. 


60 


0. 


.17 


72 


.45 


78 


.07 


2 


.464 


-28. 


.5 


31 


.0 


2 


.12 


13 48 


20 


89 


+30 


20 


05 


.3 


1958. 


.373 


1996. 


515 




17 


36 


18. 


.09 


-0 


.73 


16. 


66 


17. 


.57 


-0 


.91 


22 


.49 


23 


.53 


1 


.865 


-27. 


.7 


33 


.2 


1 


.72 


13 53 


26 


.08 


+36 


20 


19 


.6 


1950. 


.452 


1996. 


518 




16. 


30 


15 


.90 


0. 


40 


15. 


97 


15. 


63 





.34 


1 


.70 


1 


.57 





.286 


-25. 


.4 


30 


.5 


-0 


.12 


13 55 


29 


.94 


+29 


30 


58 


.9 


1958. 


.373 


1996. 


518 




18 


25 


18 


.47 


-0 


.22 


17. 


43 


17. 


33 


0. 


.10 


1 


.71 


1 


.51 





.302 


-23. 


.0 


30 


.6 





.91 


14 10 


36 


.79 


+29 


55 


50 


.7 


1958. 


.373 


1996. 


.518 




18 


61 


17 


.69 





92 


18. 


11 


17. 


76 





38 


4 


.00 


4 


.09 





.570 


-25. 


.2 


31 


.5 





.93 


14 15 


28 


.45 


+37 


06 


21 


.7 


1950 


447 


1996. 


518 


* 


18 


20 


18. 


.37 


-0 


.17 


17 


76 


17. 


45 





31 


406. 


.30 


409 


.69 


2 


.360 


-27 


9 


34 


.7 


3 


.04 


14 23 


26 


.09 


+32 


52 


20 


.6 


1950 


466 


1997. 


496 




17 


25 


16 


.67 





.58 


16. 


16 


16. 


23 


-0 


.07 


8 


.71 


8 


.16 


1 


.903 


-29. 


1 


32 


.8 





.72 


14 25 


50 


.71 


+24 


04 


03 


.5 


1950 


.299 


1997. 


496 




17 


93 


17 


.78 





.15 


17. 


13 


17. 


11 


-0 


.01 


121. 


14 


320 


.88 





.654 


-25. 


.4 


33 


.5 


2. 


.85 


14 31 


20 


.55 


+39 


52 


41 


.8 


1950. 


532 


1996. 


523 


* 


16 


47 


16 


.68 


-0 


.21 


15. 


87 


15. 


69 





.18 


207. 


.87 


210 


.01 


1 


.217 


-28. 





33 


.8 


2. 


.17 


14 37 


19 


.19 


+38 


01 


38 


.1 


1950. 


.466 


1996. 


690 




18 


24 


17 


.57 


0. 


67 


17 


75 


17. 


41 





.34 


34 


.25 


35 


.08 


1 


.043 


-26. 


.8 


32 


.9 


1 


.76 


14 37 


56 


.48 


+35 


19 


36 


.8 


1950. 


.466 


1996. 


526 




18. 


16 


17. 


.85 





31 


17. 


95 


17. 


.42 





.53 


16 


.53 


53 


.48 





.537 


-24. 


.9 


32 


.5 


2. 


.12 


14 55 


43 


.46 


+30 


03 


22 


.6 


1955. 


364 


1996. 


.282 




17. 


39 


16 


.69 


0. 


70 


16. 


94 


16. 


21 





.73 


1 


.80 


1 


.41 





.626 


-26. 


.4 


31 


.2 





.17 


15 14 


43 


.03 


+36 


50 


50 


.4 


1956 


.499 


1996. 


.282 




16 


.49 


15 


.78 





.71 


16. 


30 


15. 


87 





.43 


48 


.97 


70 


.99 





.369 


-26. 


1 


32 


.3 


1 


.14 


15 17 


28 


.49 


+28 


56 


15 


.8 


1954. 


.564 


1996. 


282 




18 


65 


17. 


.86 


0. 


79 


17. 


53 


17. 


40 





.13 


1 


.39 





.81 





.208 


-22. 


.7 


30 


.1 





.59 


15 19 


36 


.16 


+28 


38 


27 


.0 


1954 


564 


1996. 


.282 




17 


22 


17 


.37 


-0 


.15 


16. 


.61 


16. 


75 


-0 


.14 


2 


.08 


1 


.99 





269 


-23. 


,8 


30 


.5 





.56 


15 23 


11 


.50 


+37 


59 


28 


.4 


1950. 


.447 


1996. 


690 




18 


29 


18. 


.14 


-0 


.15 


17. 


43 


17. 


64 


-0 


.21 


1 


.67 


1 


.83 


1 


.344 


-26. 


.5 


31 


.9 





.80 


15 25 


23 


.58 


+42 


01 


17 


.1 


1955. 


329 


1996. 


523 




18 


43 


18 


.47 


-0 


.04 


17 


71 


17. 


61 





.10 


104. 


.19 


122 


19 


1 


.198 


-26. 


.2 


33 


6 


2. 


.65 


15 29 


18 


.04 


+32 


48 


42 


.3 


1950 


447 


1996. 


690 




17 


.70 


17 


.52 





18 


17. 


05 


16. 


99 


0. 


.06 


28 


.91 


29 


.62 


1 


652 


-27. 


9 


33 


.2 


1 


.61 


15 30 


45 


.79 


+38 


39 


52 


.7 


1955. 


329 


1996. 


282 




18 


54 


17 


.91 


0. 


63 


17. 


82 


17. 


36 





.46 


2 


.28 


2 


.28 


2 


.035 


-28. 


.0 


32 


.3 





.62 


15 37 


30 


.98 


+33 


58 


36 


.1 


1950 


447 


1996. 


696 




18 


.34 


17 


.74 





60 


17. 


22 


17. 


43 


-0 


.21 


5 


.27 


5 


.58 


1 


.025 


-26. 


.6 


32 


.1 


1 


.03 


15 40 


42 


.98 


+41 


38 


16 


.0 


1955 


.329 


1996. 


696 




17 


70 


17. 


.93 


-0 


.23 


17. 


28 


17. 


15 


0. 


13 


17 


.67 


18 


.28 


2 


.506 


-28. 


,5 


33 


.1 


1 


.51 


15 42 


27 


.03 


+29 


42 


02 


.4 


1950. 


460 


1996. 


282 




18 


55 


18. 


10 


0. 


45 


17. 


88 


17. 


61 





.24 


1 


.54 


1 


.35 





.457 


-24. 


3 


30 


.9 





.69 


15 42 


41 


.12 


+29 


31 


28 


.9 


1950. 


.460 


1996. 


282 




18. 


68 


18. 


.19 


0. 


49 


17. 


99 


17. 


90 





09 


1 


.58 


1 


.81 


1 


.098 


-26. 


3 


31 


.7 





.72 


15 43 


48 


.62 


+40 


13 


24 


.8 


1955. 


329 


1997. 


578 




17. 


91 


17. 


.72 





19 


17. 


11 


16. 


79 





.35 


3 


.24 


2 


.75 





.318 


-23. 


.8 


30 


.9 





.88 


15 44 


05 


.61 


+32 


10 


49 


.1 


1950. 


.447 


1997. 


.578 




17. 


70 


17 


.92 


-0 


.22 


17. 


01 


17. 


52 


-0 


.48 


210. 


.79 


213 


.88 


1. 


.047 


-26. 


.4 


33 


.7 


2 


.69 


15 45 


12 


.93 


+30 


05 


08 


.0 


1950 


.460 


1996. 


523 




18. 


50 


19 


.24 


-0 


.74 


17. 


.71 


17. 


83 


-0 


.12 


6 


.73 


6 


.26 





.322 


-22. 


3 


31 


.2 


1 


.81 


15 48 


17 


.91 


+35 


11 


27 


.5 


1950. 


.359 


1996. 


523 




18 


26 


18. 


.04 





22 


17. 


92 


17. 


70 





22 


140. 


.94 


141 


.51 





479 


-24. 


.5 


32 


.9 


2. 


.63 


15 54 


29 


.39 


+30 


01 


19 


.0 


1950 


.460 


1996. 


.523 


* 


18 


.38 


18. 


.08 





.30 


17. 


75 


17. 


18 





.57 


40 


.04 


41 


.22 


2 


.686 


-28. 


6 


33 


.8 


1 


.91 


15 58 


55 


.20 


+33 


23 


19 


.4 


1950 


359 


1996. 


.523 




18. 


02 


17. 


.96 





06 


16. 


94 


16. 


91 


-0 


.00 


142. 


.96 


144 


.89 


1. 


.655 


-27. 


.5 


33 


.9 


2. 


.18 


16 02 


57 


.37 


+30 


38 


51 


.6 


1950. 


.512 


1996. 


523 




18. 


05 


18. 


.24 


-0 


.19 


17. 


60 


17. 


81 


-0 


.21 


1 


.23 


1 


.39 





.809 


-25. 


.5 


31 


.3 





.65 


16 03 


51 


15 


+30 


02 


08 


.5 


1950. 


512 


1996. 


.523 




18. 


17 


17. 


.88 


0. 


29 


17. 


72 


17. 


50 





.22 


53 


.69 


54 


.18 


2 


.026 


-28. 


.1 


33 


.7 


1 


.99 


16 06 


48 


.11 


+29 


10 


48 


.0 


1950 


512 


1996. 


523 




18. 


.82 


18. 


.12 





40 


17. 


67 


17. 


40 





.27 


1 


.26 





.74 





.328 


-23. 


.2 


30 


.5 





.75 
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Long-Term Quasar Variability 
Table 2 — Continued 



RA 


Dec 


POSS1 obs 


CCD obs 


B 


O 


B-O 


R 


E 


R-E 


Sp Si 


z M(B) 


L(R) 


log R* 


Morpholo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) (12) 


(13) (14) 


(15) 


(16) 


(17) 



16 


11 


23 


26 


+29 


59 47.2 


1950 


512 


1996 


526 


18 


50 


18 


56 


-0 


06 


17 


85 


17 


75 





10 


1 


69 


1 


90 


1 


711 


-27 





32 


1 





83 


16 


13 


41 


09 


+34 


12 48.3 


1954 


564 


1996 


526 


18 


11 


17 


74 





37 


17 


23 


17 


53 


-0 


30 


3532 


Od 


3605 


54 


1 


404 


-27 


3 


35 


2 


3 


81 


16 


14 


36 


84 


+28 


39 06.2 


1950 


512 


1996 


523 


19 


39 


19 


05 





3d 


17 


80 


17 


51 





29 


2 


22 


1 


31 





316 


-22 


5 


30 


7 


1 


25 


16 


19 


02 


43 


+30 


30 51.4 


1950 


512 


1996 


526 


17 


24 


17 


29 


-0 


05 


16 


56 


16 


46 





10 


66 


82 


67 


58 


1 


284 


-27 


5 


33 


d 


1 


91 


16 


23 


07 


39 


+30 


04 06.2 


1954 


567 


1996 


526 


18 


08 


17 


29 





79 


17 


32 


16 


95 





37 


1 


65 


1 


Od 


1 


165 


-27 


3 


31 


7 





31 


16 


23 


19 


92 


+41 


17 02.8 


1953 


460 


1997 


496 


17 


d6 


17 


32 





14 


16 


64 


16 


67 


-0 


03 


3 


23 


2 


98 


1 


618 


-28 


1 


32 


3 





57 


16 


23 


30 


56 


+35 


59 32.8 


1954 


564 


1996 


526 


18 


55 


17 


91 





64 


17 


90 


17 


59 





31 


259 


58 


266 


83 





867 


-26 





33 


6 


2 


80 


16 


24 


22 


02 


+39 


24 41.5 


1953 


460 


1996 


526 


18 


94 


18 


31 





63 


18 


26 


17 


76 





50 


132 


32 


144 


52 


1 


117 


-26 


2 


33 


6 


2 


66 


16 


25 


48 


78 


+26 


46 58.5 


1954 


567 


1997 


595 


17 


40 


17 


09 





31 


17 


01 


17 


05 


-0 


Od 


10 


12 


9 


69 


2 


521 


-29 


4 


33 


1 





91 


16 


26 


59 


26 


+30 


15 34.9 


1954 


567 


1996 


526 


18 


36 


18 


17 





19 


17 


55 


17 


19 





36 


5 


34 


5 


42 


1 


582 


-27 


2 


32 


5 


1 


14 


16 


29 


01 


32 


+40 


08 00.0 


1953 


d60 


1996 


526 


17 


67 


18 


11 


-0 


dd 


17 


25 


17 


86 


-0 


61 


11 


97 


11 


94 





271 


-23 


1 


31 


3 


1 


61 


16 


30 


20 


76 


+37 


56 56.2 


1954 


564 


1996 


690 


16 


59 


17 


31 


-0 


72 


16 


43 


17 


45 


-1 


02 


20 


00 


21 


73 





395 


-24 


7 


31 


9 


1 


53 


16 


33 


02 


11 


+39 


24 27.5 


1953 


460 


1996 


690 


17 


17 


16 


86 





31 


16 


56 


15 


97 





59 


41 


29 


54 


38 


1 


021 


-27 


4 


33 


1 


1 


67 


16 


33 


48 


99 


+31 


3d 11.8 


1954 


567 


1996 


690 


17 


38 


17 


72 


-0 


3d 


16 


87 


17 


27 


-0 


40 


1 


79 


2 


07 


1 


516 


-27 


5 


32 








55 


16 


34 


02 


97 


+39 


00 00.7 


1953 


460 


1996 


690 


18 


91 


18 


38 





53 


18 


22 


18 


03 





19 


915 


40 
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Note. — Descriptions of table columns: 
Cols 1-2: FIRST radio positions (J2000). 
Col 3: Date of POSS-I observation. 
Col 4: Date of CCD observation. 

Cols 5-10: Extinction-corrected CCD magnitudes (B, R), POSS-I magnitudes (O, E) and colors. An asterisk preceding B indicates sources with too 
CCD calibration stars to recalibrate the POSS-I magnitudes (see text.) 

Cols 11-12: 20 cm peak (col. 11) and integrated (col. 12) flux densities in mjy from the FIRST catalog. 

Col 13: Redshift. 

Col 14: Absolute B magnitude. 

Col 15: Log of radio luminosity (erg cm~ 2 s _1 Hz -1 ) at a rest-frame frequency of 5 GHz, assuming spectral index a = —0.5 (F v oc v a ). 
Col 16: i<"-corrected 5 GHz radio to 2500 A optical luminosity ratio using the definition of Stocke et al. (1992). 

Col 17: Radio morphology. Categories are FRII (double, with approximate separation), core-jet (with approximate jet length), and complex. Objects w 
no notation are symmetrical, isolated sources. 



